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ABSTRACT 

Context. Although oxygen is an important tracer of Galactic chemical evolution, measurements of its abundance in the atmospheres 
of the oldest Galactic stars are still scarce and rather imprecise. This is mainly because only a few spectral lines are available for the 
abundance diagnostics. At the lowest end of the metallicity scale, oxygen can only be measured in giant stars and in most of cases 
such measurements rely on a single forbidden [O i] 630 nm line that is very weak and frequently blended with telluric lines. Although 
molecular OH lines located in the ultraviolet and infrared could also be used for the diagnostics, oxygen abundances obtained from 
the OH lines and the [Oi] 630 nm line are usually discrepant to a level of ~ 0.3 - 0.4 dex. 

Aims. We study the influence of convection on the formation of the infrared (IR) OH lines and the forbidden [Oi] 630 nm line 
in the atmospheres of extremely metal-poor (EMP) red giant stars. Our ultimate goal is to clarify whether a realistic treatment of 
convection with state-of-the-art 3D hydrodynamical model atmospheres may help to bring the oxygen abundances obtained using the 
two indicators into closer agreement. 

Methods. We used high-resolution (R = 50 000) and high signal-to-noise ratio (S /N ~ 200 - 600) spectra of four EMP red giant 
stars obtained with the VLT CRIRES spectrograph. For each EMP star, 4-14 IR OH vibrational-rotational lines located in the spectral 
range of 1514-1548 and 1595-1632 nm were used to determine oxygen abundances by employing standard ID local thermodynamic 
equilibrium (LTE) abundance analysis methodology. We then corrected the ID LTE abundances obtained from each individual OH 
line for the 3D hydrodynamical effects, which was done by applying 3D-1D LTE abundance corrections that were determined using 
3D hydrodynamical CO^BOLD and ID hydrostatic LHD model atmospheres. 

Results. We find that the influence of convection on the formation of [O i] 630 nm line in the atmospheres of EMP giants studied here 
is minor, which leads to very small 3D-1D abundance corrections (Asd-id ^ -0.01 dex). On the contrary, IR OH lines are strongly 

affected by convection and thus the abundance corrections for these lines are significant, A 313 - 1 D ~ -0.2-0.3 dex. These abundance 

corrections do indeed bring the ID LTE oxygen abundances of EMP red giants obtained using IR OH lines into better agreement with 
those determined from the [O i] 630 nm line. Since in the EMP red giants IR OH lines are typically at least a factor of two stronger 
than the [O i] line, OH lines may be useful indicators of oxygen abundances in the EMP stars, provided that the analysis is based on 
3D hydrodynamical model atmospheres. 

Key words. Stars: Population II - Stars: late type - Stars: atmospheres - Stars: abundances - Techniques: spectroscopic - Convection 
- Hydrodynamics 


1. Introduction 

Oxygen is the most abundant element in the Universe, after hy¬ 
drogen and helium, a nd is the main product of the Type II SNe 
(see e.g. lArnetlll996l. and references therein). Massive stars have 
lifetimes of the order of 10® to 10^ years and end their lives as 
sup ernovae, releasin g oxygen into the interstellar medium (see 
e.g. lMatteuccill201^ and references therein). As a first approx- 


* Based on observations obtained at the European Southern Observa¬ 
tory (ESO) Very Large Telescope (VLT) at Paranal Observatory, Chile 
(observing programme 089.D-0079). 


imation it may be expected that the oxygen abundance in ex¬ 
tremely metal-poor (EMP) stars correlates with their age, the 
most oxygen-poor stars being the oldest. On closer consideration 
it becomes clear that the oxygen abundance at any given time 
also depends on the star formation rate. In a complex system, 
like our Galaxy, that may host populations of different origins 
(perhaps some accreted by dwarf satellite galaxies) the evolution 
of the O/Ee ratio with time depends on the details of the chemi¬ 
cal evolution of the different populations. The abundance ratio of 
oxygen to other elements, and in particular to iron, places strong 
constraints on proposed scenarios of Galactic chemical evolu- 
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tion. The precise knowledge of oxygen abundances in the oldest 
Galactic populations is important since it allows us to discrimi¬ 
nate among different models (see lMatteuc(^l2012L chapter 5.3, 
for a discussion). 

In spite of these virtues, the oxygen abundance is actu¬ 
ally very difficult to measure in the metal-poor stellar atmo¬ 
spheres; the only indicators available in the EMP stars are (i) 
the forbidden [Oi] line at 630 nm for giants, (ii) the permit¬ 
ted Oi IR triplet for dwarfs (777.2 - 777.5 nm), and (iii) OH 
lines in the UV and IR, for both giants and dwarfs. Unfortu¬ 
nately, these different indicators often provide discrepant oxy¬ 
gen abundances, which may be due to non-local thermody¬ 
namic equilibrium (NLTE) effects (especially for the IR triplet), 
3D hydrodynamical effects (especially for the OH lines but 
also for the triplet), incorrectly determined atmospheric parame- 
ters (e.g. JOO is very sensitive to the surfa ce gravity; see e.g. 
ICavrel etdl 1200 It iGarcfa Perez et al.l l2006b , and so on. Con¬ 
sequently, reliable oxygen abundances of the EMP stars are 
scarce: the number of objects in which they have been mea¬ 
sured so far is small while the measurements themselves are fre¬ 
quently not sufficiently reliable. This has led different groups to 
claim the existence of different trends of [O/Eel with decreasing 
metallicity (see e.g. Boesgaard et al.lfTPOQt llsraelian et al.ll200ll 


lAkerman et al.ll2064ll . which makes discriminating between the 

various proposed scenarios of Galactic oxygen evolution very 
difficult. Clearly, such a situation is unsatisfactory and calls for 
the re-assessment of oxygen abundances in EMP stars. 

It should be noted that oxygen abundances in the EMP gi¬ 
ants are nearly always obtained from the measurements of a 
single forbidden [Oi] 630 nm line. One reason why OH lines 
have rarely been used for the abundance determinations so far 
( Melendez et alJl200ll : iMelendez & Barbu^l2002l : iBarbuv et alJ 
I 2 OO 3 I 1 is that the oxygen abundances obtained from these lines 
are signi ficantly different fro m those determined using the [O i] 
line (e.g. lBarhuv et al.l[2003h . Part of the problem may be related 
to the [Oi] line measurements, because this line is very weak 
and often blended with telluric lines. On the other hand, molec¬ 
ular lines may be prone to 3D hydrodynamical effects, such as 
horizontal temperature fluctuations, which may lead to substan¬ 
tial differences in the oxygen abundances derived with the 3D 
hydrodynamical and classical ID hydrosta tic models (see e.g. 
ICollet et aDl2007l ; lDobroyolskas et al.ll2013h . It may thus be rea¬ 
sonable to anticipate that oxygen abundances obtained from OH 
lines and those determined using the [O i] line could be brought 
closer if these effects are properly taken into account in the 
model atmospheres used in the abundance determination. If this 
were proyed to be true, OH lines could be used as independent 
yaluable indicator of the oxygen abundance, with the adyantage 
that (i) OH lines are numerous in the infrared and thus one would 
not need to rely on a single line for measuring the oxygen abun¬ 
dance, as in the case with the forbidden [O i] 630 nm line; and (ii) 
infrared (IR) OH lines are at least a factor of 2 stronger than the 
[O i] 630 nm line is, which makes them easier to measure and 
may lead to a better accuracy in the determined oxygen abun¬ 
dances. 

In this pilot study we therefore inyestigate whether such con¬ 
siderations are indeed yalid in the case of EMP red giant stars. 
Eor this purpose, we used high-resolution and high signal-to- 
noise spectra of four EMP giants obtained with the CRIRES 
spectrograph at the VET UTl telescope to determine ID LTE 
oxygen abundanc es using IR OH lin es. Then, 3D hydrody- 
namical CO^BOLD (jEreytag et alj|201^ and ID hydrostatic LHD 
(ICaffau & Ludwig i2007h model atmospheres were used in order 
to assess the role and influence of conyection on the formation 


of IR OH lines and the forbidden [O i] line, as well as to deter¬ 
mine the 3D-corrected abundances of oxygen from IR OH lines. 
Einally, we discuss the possible role of yarious other effects that 
may influence the formation of IR OH lines and that may need 
to be taken into account in the abundance determinations using 
IR OH lines. 


2. Methodology 

2.1. Observations and data reduction 

Since IR OH lines are yery weak in EMP giants, one needs to 
haye access to spectra with yery high S /N ratios in order to ob¬ 
tain reliable oxygen abundances. Eor this reason, we limited our 
target list to objects with V < 12 . 0; the se are the brightest EMP 
giants in the list of ICayrel et alJ (l2004l) . This has allowed us to 
keep the obserying time request within reasonable limits eyen 
in the case of the faintest star, CS 31082-001; we note that ex¬ 
cept for HD 122563, no EMP giants haye been obseryed with 
CRIRES so far. The list of our EMP targets is proyided in Ta¬ 
ble [U 

High-resolution //-band near-infrared spectra of the four 
EMP giants were obtained in s eryice mode with the CRIRES 
spectrograph dKaufl et alJl2004ll located at the VET UTl tele¬ 
scope (programme ID 089.D-0079(A), PI; A. Kucinskas). We 
used the entrance slit width of 0.4 arcsec for our obseryations 
and did not use the adaptiye optics system. The spectra were 
acquired in ABBA nodding pattern for dark current and sky 
emission remoyal, using two CRIRES settings: order 37 with 
the central wayelength of 1536.4 nm and a spectral range of 

1513.4 - 1548.6 nm; and order 35 with the central wayelength of 

1619.4 nm and a spectral range of 1594.9 - 1632.5 nm. The tar¬ 
gets were obseryed with a spectral resolution of /? = 50 000 and 
signal-to-noise ratio per pixel S /N ^ 200 - 600. Standard steps 
inyolyed in the reduction of the raw spectra (bad pixel correc¬ 
tion, detector non-linearity correction, flat-fielding, and waye¬ 
length calibration) were performed using the CRIRES pipeline 
(y.2.3.1) in the command-line (EsoRex) mode. 

Unfortunately, the wayelength calibration solution proyided 
by the CRIRES pipeline had systematic shifts which made spec¬ 
tral line identification difficult. Therefore, we performed a cus¬ 
tom wayelength calibration using telluric lines present in the 
obseryed spectral range. Eor this purpose, we computed the 
synthetic telluric sp ectrum using the TAPAS online seryice Q 
(iBertaux et al.l 120141) for the conditions matching those of our 
obseryations (location, date, time, and air-mass). We did not at¬ 
tempt to clean stellar spectra from the telluric lines, thus for fur¬ 
ther spectroscopic analysis we used only those OH lines which 
were free from telluric contamination. In fact, we obseryed a tel¬ 
luric star HD 24388, but because of a significantly lower S/N 
ratio this spectrum was not used to correct the target star spectra 
for the telluric features. Instead, we chose the synthetic telluric 
spectrum computed using TAPAS to identify telluric lines in our 
CRIRES spectra and to ayoid them in further abundance analy¬ 
sis. Continuum normalization of the o bseryed CRIRES sp ectra 
was made using the DECH20T package (iGalazutdinovI 1 99^ . 

2.2. Modei atmospheres 

Three types of model atmospheres were employed in this study. 
Eirst, ID LTE oxygen abundances were determined b y using 
model atmospheres computed with the ATLAS 9 code (iKuruczl 

* http://ether.ipsl.jussieu.fr/tapas/ 
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Table 1. List of EMP giants observed with CRIRES. 


Star 

a(J2000) 

b(J2000) 

V, 

(mag) 

H, 

(mag) 

^exp? 

(s) 

S/N 

Teff, 

(K) 

log0 

[cgs] 

(kms-') 

[Fe/H] 

HD 122563 

14 02 31.85 

-H094109.9 

6.2 

3.8 

35 

560 

4598 

1.60 

2.0 

-2.82 

HD 186478 

19 45 14.14 

-17 29 27.1 

9.2 

6.6 

270 

440 

4700 

1.30 

2.0 

-2.59 

BD-18:5550 

19 58 49.74 

-18 1211.1 

9.4 

6.7 

270 

600 

4750 

1.40 

1.8 

-3.06 

CS 31082-001 

01 2931.13 

-16 0045.4 

11.3 

9.6 

1500 

190 

4825 

1.50 

1.8 

-2.91 


Note; at mospheric parameters and iron abundances of the target stars a re from Spite et ahl (l2005h . except for HD 122563 which 
are from ICreevev et af] (l2012h . Microturbulence velocities, are from lSnite et alJ ( 2005 ). V and H magnitudes are from SIM- 
BAD. 


Table 2. Oxygen abundances in EMP giants. 


Star 

A(0)iumh,a 

OH 

o-A(0)“^'^"''" 

OH 

A(0)-^'JLTb,a 

OH 

crA(0)^'^'^"^'“ 

OH 

A(0)l^l^lLb 

[Oi] 

[Oi] 

HD 122563 

6.63 

0.10 

6.39 

0.11 

6.57 

0.15 

HD 186478 

7.11 

0.09 

6.85 

0.08 

6.93 

0.14 

BD-18:5550 

6.73 

0.13 

6.48 

0.15 

6.13 

0.26 

CS 31082-001 

6.98 

0.13 

6.73 

0.10 

6.46 

0.20 


“ determined in this work; from ICavrel et al.l (l2004h 


Il993h. We used th e ATLA S9 v ersion that was po rted to Linux 
by Isbordone et all (l2004h and ISbordon^ (l2005h . The models 
were calculated using opa city distribution functions (ODFNEW) 
from ICastelli et al.l (l2003ll : the latter were computed with the mi¬ 
croturbulence velocity of ^ = 2 km s '. All ATLAS9 model at¬ 
mospheres were calculated using the mixing length parameter 
Q'mlt = 1.25, with overshooting switched off. 

In addition, 3D hydrodynamical and ID hydrostatic model 
atmospheres were used for computing 3D-1D abundance cor¬ 
rections (see Sect. \32\ . For this purpose we employed a sin¬ 
gle 3D hydrodynamical model atmos phere computed usin g the 
CO^BOLD model atmosphere package (iFrevtag et al.ll201^ with 
the following atmospheric parameters: Tefr = 4595 K, loggr = 
1.6, [M/H] = -2.5 (Klevas et al., in prep.). Our choice was de¬ 
termined by the fact that 3D hydrodynamical models are rather 
expensive to compute and thus they are not yet available for 
all combinations of atmospheric parameters. Nevertheless, atmo¬ 
spheric parameters of this particular model atmosphere are suf¬ 
ficiently similar to those of FMP giants studied here: even in the 
most extreme cases the differences do not exceed ATefr = 230 K, 
Aloggf = 0.3, and A [M/H] = 0.5, while generally they are even 
smaller (see Table[T]i. 


The 3D hydrodynamical model atmosphere was computed in 
Cartesian geometry with a spatial resolution of 220 x 220 x 280 
cells in the x,y,z directions, respectively, which corresponds 
to spatial dimensions of 3.85 x 3.85 x 2.21 Gm. The hori¬ 
zontal size of the model was large enough to accommodate 
» 10 granules which allowed convective d ynamics to be reli¬ 
ably captured in the s tellar atmosphere (see Collet et ^ l2007[ 


iLudwig & KucinskasI 1201 2[ iKucinskas et al 


l2013h . We used 


opacities from the MARCS model atmosphere package, which - 
in order to make the model computations faste r - were grouped 
into six opaci t y bins using opacity binning (INordlundl [l982t 
iLudwigI Il992t Ludwig et S 1994: Vogler et al.| boO^ Solar 
chemical composition from Grevesse & Sauvall (11998 ) with a- 
element enhancement of -1-0.4 dex was used in the model calcu¬ 
lations. For CNO elements, we used solar values of A(C) = 8.41 , 
A(N) = 7.8, and A(0) = 8.66 taken from lAsnlund et all (l2005h . 


which were then scaled down to the metallicity of our 3D hy¬ 
drodynamical model atmosphere with a-element enhancement 
of h- 0.4 dex. 

In order to speed up the 3D spectral line synthesis calcula¬ 
tions, we selected 20 3D model structures computed at different 
instants in time (snapshots). When making the snapshot selec¬ 
tion we tried to ensure that the average properties of the selected 
snapshot subsample would remain similar to the average proper¬ 
ties o f the full 3D model simulation run (see e.g. IKucinskas et akl 
1201 3l for the details on the snapshot selection procedure). This 
snapshot subsample was used in subsequent 3D spectral synthe¬ 
sis computations. 

The 20-snapshot ensemble was also used to compute aver¬ 
age (3D) model, which was obtained by averaging the 3D struc¬ 
tures in the 20-snapshot subsample on surfaces of equal optical 
depth. As has been already discussed in our previous studies, 
such average (3D) models are useful in evaluating the role of 
horizontal temperature fluctuations in the spectral line forma¬ 
tion. Since the average (3D) models do not contain information 
about the horizontal fluctuations of thermodynamic properties, 
the 3D-(3D) differences provide important information about 
the role of these fluctuations in the spectral line for mation (see 
e.g. lDobrovolskas et akll^oT^ IKucinskas et aO2013ll . 

The ID hydrostatic model atmospheres that we used in the 
evaluation of 3D-1D abundance cor rections were com puted us¬ 
ing LHD model atmosphere package (ICaffau et al.ll2008l) . It is im¬ 
portant to note that 3D hydrodynamical and ID hydrostatic LHD 
model atmospheres employed in the computation of abundance 
corrections shared identical atmospheric parameters, chemical 
composition, opacities, and equation of state. This allowed us 
to make a strictly differential comparison and to estimate the 
impact of convection on the formation of spectral lines. There¬ 
fore, we used the (differential) 3D-1D abundance corrections 
obtained in this way to correct the abundances determined us¬ 
ing the ATLAS9 hydrostatic model atmosphere computed with 
different opacities. The ID LHD models were calculated with the 
mixing-length parameter aMLX = 1 0. 
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We note that effects of sphericity may in fact affect the 
structure of the model atmosphere, both in ID (see e.g. 
[Heiter & Erikssonl2006t) and, possibly, in 3D as well. To assess 
the size of these effects in ID, we performed a comparison of 
a spherically symmetric and a plane-parallel model atmospheres 
with Teff = 4500 K, log g = 1.5, [M/H] = -2.5. The choice 
was dictated by the available models, and in fact the spheri¬ 
cal model was a MARCS model, the plane-parallel an ATLAS9 
model. The spherical model had a radius of 29 Rq and an atmo¬ 
spheric extension of 0.5 Rq. We synthesized a few typical OH 
lines and the [Oi] line with LinforSD, under the assumption 
of Cartesian geome t ry. De spite this approximation, according to 
iHeiter & Eriksson! (l2006l) . the main effect should be captured 
by this procedure. We found a slight strengthening of the [O i] 
line (r; 0.01 dex), and a modest strengthening of the OH lines 
(ai 0.06 dex) in the spherical versus the plane-parallel model. 
From this we conclude that sphericity effects are indeed present 
but on a level significantly smaller than the 3D effects. 

Unfortunately, we are currently not able to assess the size 
of sphericity effects on the 3D hydrodynamical model struc¬ 
tures because the computation of spherical 3D hydrodynamical 
CO^BOLD models is not yet feasible. Nevertheless, we do not ex¬ 
pect that sphericity effects play a major role for the 3D-1D abun¬ 
dance corrections discussed here since the modest extension of 
the atmosphere is unlikely to change the granular dynamics no¬ 
ticeably. The overall effect of the geometry should largely cancel 
out because of the differential nature of the 3D-1D corrections. 

3. Results 

3.1. ID LTE oxygen abundances 

To determine oxygen abundances in the target EMP giants, we 
used infrared OH vibrational-rotational lines (X^H) from the 
first-overtone sequence (Ay = 2) located between 1514 - 1548 
and 1595 - 1632 nm (see Table I aTI and Figs. IB.lllB.5l in On¬ 
line Material). We note that the spin-orbit interaction gives rise 
to two spin-split substates: ^ni /2 and ^n 3 / 2 . All the lines that 
we consider are from the P-branch (A7 = -1), the rotational 
levels are subject to A-doubling, i.e. splitting that is due to per- 
turbation by the rotati onal levels of other electronic states (see 
ICoxon & Foster!Il982h . The equivalent widths of the observed 
lines were measured using the splot task in IRAF0 and are 
provided in Table lA.ll Although some of the OH lines used in 
the abundance determination are contaminated by telluric lines, 
these telluric lines are in fact very weak. Therefore, we did not 
correct for their influence when measuring the equivalent widths 
of OH lines. We note that all observed targets have known [Oi] 
line strengths, as well as oxygen abundances determined using 
this line (provided in Table |2ll. As expected, in most cases the 
IR OH lines are significantly stronger than the forbidden [O i] 
line at 630 nm. 

We also note that although there are CN lines in the spec¬ 
tral range investigated in this study, they are very weak (< 
0.05% depth from continuum) and their influence on the oxy¬ 
gen abundance determination is negligible, even in cases where 
they blend with OH lines. All the stars investigated in this 
study have a normal or low car bon abundance of [C/Fe] < 0.2 
dTakeda & Takada-Hidaill2013h . which leads to very weak CN 
lines in their IR spectra. 

^ IRAF is distributed by the National Optical Astronomy Observato¬ 
ries, which are operated by the Association of Universities for Research 
in Astronomy, Inc., under cooperative agreement with the National Sci¬ 
ence Foundation. 


Oxygen ID FTE abundances were determined using the 
curve-of-growth method, employing ATLAS 9 model atmo¬ 
spheres computed as described in Sect. 12.21 above. The ATLAS9 
models were calculated using the atmospheric parameters of the 
individual EMP giants that were taken from the literature and 
are provided in Table [1] Synthetic spectra used in the abun¬ 
dance determination were computed usin g the Finux version of 
the SYNTHE spe ctral synthesis package (ISbordone et al.l 120041 : 
ISbordonel2005l) . We then measured the equivalent widths of the 
OH lines in the synthetic spectra using the splot task in IRAF, 
i.e. using the same methodology as for the observed lines. Fi¬ 
nally, ID FTE oxygen abundances were determined by interpo¬ 
lating between the curves of growth at the equivalent widths of 
IR OH lines that we measured in each EMP giant (Table|2|. 

We estimate that the uncertainty in the determined oxygen 
abundance arising due to the uncertainty in the continuum level 
determination (which we estimate to be +0.2 %) is ±0.06 dex. 
The uncertainty in the measured oxygen abundance provided in 
Table |2] is a square root of the sum of (i) the dispersion due to 
line-to-line variations of the determined oxygen abundance for 
each EMP star and (ii) the squared uncertainty due to continuum 
placement. 

Many OH lines studied here form from transitions that in¬ 
volve rotational levels which, because of perturbations by other 
electronic states, are split into two sub-levels (A-doubling). 
Therefore, the two lines arising because of such transitions 
should give identical oxygen abundances. Since in reality this, 
apparently, is not the case, the difference between abundances 
obtained from the two individual components may provide ad¬ 
ditional constraints on the size of the uncertainties in the deter¬ 
mined oxygen abundances. In our case, we find that the largest 
difference in the oxygen abundance obtained from two lambda¬ 
doubling components is ~ 0.12 dex, which is comparable with 
the uncertainties provided in Table 

3.2. 3D-1D abundance corrections and 3D LTE oxygen 
abundances 

The 3D FTE abundances of oxygen were obtained by using 3D- 
ID abundance corrections. The 3D-1D abundance correction, 
A 3 D-ID = A(Xi) 3 D - A(Xi)iD, is the difference in the abun¬ 
dance of element Xi that would be obtained at a fixed equiva¬ 
lent width of a given spectral line with 3D hydrodynamical and 
ID hydrostatic model atmospheres, respectively. In this work, 
abundance corrections were obtained using the curves of growth 
computed with 3D hydrodynamical CO^BOLD and ID hydrostatic 
LHD model atmospheres. Synthetic spectral line profiles were 
calculated under the assurnption of FTE in 3D and ID, in both 
cases using the LinforSlQ spectral synthesis package. Abun¬ 
dance corrections were then computed by measuring the differ¬ 
ence between the curves of growth obtained with the 3D hydro¬ 
dynamical and ID hydrostatic model atmospheres, respectively. 

The 3D-1D abundance corrections were determined for all 
IR OH lines used in this work, as well as for the forbidden [O i] 
630 nm line, using the line parameters provided in Table lA.ll 
and Tables lBTlIB.TI Apart from A 3 D- 1 D correction, we also com¬ 
puted two additional corrections, A 3 d-< 3 d> = A(Xi) 3 D-A(Xi)( 3 D), 
and A( 3 d)-id = A(Xi)( 3 D) - A(Xi)iD. The former is useful in 
assessing the importance of horizontal fluctuations of thermo¬ 
dynamic and hydrodynamical quantities in 3D hydrodynamical 
model atmospheres, while the latter allows the role of differences 
to be estimated between the temperature profiles in the average 

^ http://www.aip.de/~mst/Linfor3D/linfor_3D_manual.pdf 
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Fig. 1. 3D LTE oxygen abundance in HD 122563, HD 186478, BD - 
18:5550, and CS 31082-001 (top to bottom), as determined from IR OH 
lines, shown as the function of reduced equivalent width. 


Fig. 2. 3D LTE oxygen abundance in HD 122563, HD 186478, BD - 
18:5550, and CS 31082-001 (top to bottom), as determined from IR OH 
lines, shown as the function of excitation potential. 


(3D) and ID model atmospheres (see Sect. |4ll. Abundance cor¬ 
rections for IR OH lines used in this work with all EMP giants 
are provided in Tables IB.lljB.4l In these tables we also provide 
abundance corrections for the forbidden [O i] 630 nm line. 

The 3D LTE oxygen abundances were obtained by adding 
A 3 D- 1 D correction to ID LTE oxygen abundances obtained in 
Sect. 13. II above. Although we are clearly aware of the many lim¬ 
itations of such approach (see Sect. |4] for discussion), we nev¬ 
ertheless think that the 3D-corrected oxygen abundances, when 
compared to those derived in ID LTE, may provide important 
qualitative information for the assessment of the role of con¬ 
vection in the spectral line formation taking place in the atmo¬ 
spheres of EMP giants. The obtained 3D LTE abundances are 
listed in Table |2] for all EMP giants studied here. The obtained 
abundances are plotted in Fig. [T] versus the line reduced equiv¬ 
alent width, logiEW!A), and in Fig. versus the line excitation 
potential, x, for all four EMP stars studied here. 

We point out that our 3D-1D oxygen abundance correc¬ 
tions for different programme stars were computed using a single 
3D hydrodynamical model atmosphere, whereas the atmospheric 
parameters of the programme stars differ by up to +230 K 
and -0.3 dex in Teff and log g with respect to those of the 3D 
model atmosphere. With the ID ATLAS9 models, this would 
amount to differences in the determined oxygen abundance of 
up to +0.65 dex and +0.12 dex, r espectively. These nu mbers are 
very similar to those obtained bv lBarbuv et alJ (l2003h . Unfortu¬ 
nately, our 3D hydrodynamical model grid is too sparse to esti¬ 


mate the size of temperature and gravity sensitivity of the 3D- 
ID abundance corrections. Nevertheless, our crude estimate is 
that an increase in the effective temperature by +230 K (from 
Tes - 4600 K) would change the full A 3 d_id abundance cor¬ 
rection by -0.06-0.08 dex (depending on the particular OH 

line). 

The OH lines investigated in this work are weak and situated 
on the linear part of the curve of growth. It makes the spectral 
lines under consideration insensitive to the choice of the micro¬ 
turbulence velocity value. Changing the microturbulence veloc¬ 
ity by 1 km s ' from the values listed in Table [T] would change 
3D-1D abundance correction and the average oxygen abundance 
by less than 0.01 dex. 

4. Discussion 

Abundance corrections provided in Tables lBTT4B.4l clearlv show 
that horizontal temperature fluctuations caused by convective 
motions in the stellar atmosphere have significantly larger im¬ 
pact on the strengths of IR OH lines than the influence of the 
different temperature stratifications in the average (3D) and ID 
model atmospheres. This is indicated by the different size of 
A 3 d-< 3 d> and A< 3 d>-id abundance corrections; the former is al¬ 
ways larger than the latter (i.e. in absolute value). The phys¬ 
ical reasons behind such behaviour are in fact rather obvious. 
As can be seen in Fig. [3 horizontal temperature fluctuations 
(as measured by ATrms = where de- 
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Fig. 3. Top: Temperature stratifications in the 3D hydrodynamical (grey 
scales indicating the temperature prohahility density), average (3D) 
(solid red line), and ID LHD (dashed red line) model atmospheres com¬ 
puted with Teff = 4600 K, log g = 1.6, [M/H] = -2.5. Horizontal 
bars show the optical depth intervals where 90% of the OH line (A = 
1605.5 nm,;i' = 0.48 eV) and [Oi] line (A = 630.0 nm,^ = 0.00 eV) 
equivalent width is formed (i.e. 5% to 95%): black and dashed blue 
bars correspond to the line forming regions in the full 3D and ID model 
atmospheres (the equivalent widths of OH and [Oi] lines are 1.6 pm 
and 0.4 pm). Bottom: RMS value of horizontal temperature fluctua¬ 
tions, AT^rms, in the 3D model (black line) and temperature difference 
between the average (3D) and ID models (dashed blue line). 



notes temporal and horizontal averaging on surfaces of equal 
optical depth, while To = {T)x,y,t, is the depth-dependent av¬ 
erage temperature) become increasingly larger in the outer at¬ 
mosphere, reaching ATrivis ~ 800 K at logrRoss ~ -5. Since 
IR OH lines form rather high in the atmosphere, this leads to 
large non-linear horizontal variations of OH number densities 
and, thus, line opacities in the line forming regions. The result is 
that IR OH lines become stronger in 3D, which in turn leads to 
negative (and relatively large) A 3 d~< 3 d> abundance corrections. 
On the other hand, the temperature difference between (3D) and 
ID is small (and mostly positive) at all optical depths, which 
results in weaker lines in (3D) and, thus, small (and positive) 
A{3d>-id corrections. 

In contrast, the influence of convection on the formation of 
the [Or] oxygen line at 630 nm is very small. This is in line 
with our earlier findings (see e.g. lDobrovolskas et alJl201^ and 
is quite easily understood. Neutral oxygen is a majority species 
in the regions where the [O i] 630 nm line forms, and because of 
its high ionization potential its abundance is insensitive to tem¬ 
perature fluctuations. At the same time, differences between the 
T profiles in the average (3D) and ID models in the line for¬ 
mation regions are small (Fig.O. This results in small A 3 d-( 3 d> 
and A( 3 d>-id abundance corrections, and thus, in small total cor¬ 
rection, A 3 D- 1 D- It is important to mention in this respect that 
the [Oil line at 630 nm is equally insensitive to NLTE effects 
(e.g lAsplund 1120051 and references therein). One may therefore 
make a relatively safe assumption that ID LTE abundances ob¬ 
tained using the [O i] 630 nm line are essentially unaffected by 



[Fe/H] 


Fig. 4. Oxygen-to-iron ratios in Galactic EMP stars plotted versus 
[Fe/H] ratio. Abundances of EMP giants studied here are marked 
by triangles: open red triangles - oxygen abundances in the pro¬ 
gramme stars determined determined in this work with classical ID 
model atmospheres; filled red triangles - oxygen abundances cor¬ 
rected for 3D hydrodynamical effects. Solid grey circles show oxy¬ 
gen abundances in red giants, subgiants, and mai n sequence stars de¬ 
termi n ed using the forbidden [Oi ] 630 n m line (|Gratton_&T)rtolani 


198f: iBarbuv & E r delvi-Mendes I Il989l; faaf tetal. 1992|j_teaAuv 

tr n • 


19881: ISneden et al.ll 19911: iFulbright^fGaflll 19991 ; ICavrel et aLll20^ 
Cowan et alJl2002l : ICavrel et alll2004h . Black triangles hi ghlight oxy- 


- . gh_ - 

:en ab undances determined from [O i] 630 nm line in ICavrel et all 
20041) for the same stars as studied here. T he solid line is the e volu- 
tionary model of the oxygen ab undance fromlFrancois et all ) l2004t) . We 
note that [O i] abundances fromICavrel et all ll2004h and the theoretical 
model from [Francois et akl l l2004h were shifted to the reference scale 
used in our study where solar oxygen abundance is A(0) = 8.66. 


the 3D and NLTE effects, and may thus serve as an important 
reference point for validating oxygen abundances obtained from 
other lines. In this respect, we find it reassuring that in our EMP 
giants the 3D LTE oxygen abundances obtained from IR OH 
lines agree well with the abundances determined from the for¬ 
bidden line (Table 121 Eig.|4|). 

We also find that measuring the IR OH lines in metal-poor 
giants is feasible at a resolution of 50 000, provided a high S/N 
ratio of a few hundred can be obtained. The lines are suffi¬ 
ciently numerous so that a few of them are always likely to be 
free from contamination of telluric lines. In this respect the syn- 
thetic spectra of the terrestrial atmosphere computed by TAPAS 
(iBertaux et al.ll20l'4l) are precious for identifying the contami¬ 
nated lines. In the present study we had enough clean lines and 
thus we felt it was not necessary to attempt to subtract the telluric 
components. Nevertheless, the fidelity of the TAPAS spectra is 
high enough so that subtraction of the telluric lines may be en¬ 
visaged in cases in which a use of a larger number of lines is 
desirable. 

Although our modelling of the formation of these lines is 
not entirely satisfactory, as can be deduced from the existence 
of slight trends of abundance with equivalent width and excita¬ 
tion potential shown in Pigs.[T] and |2] respectively, we believe 
that the use of 3D model-atmospheres improves the modelling. 
We arrive at this conclusion by noting that on average ID LTE 
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oxygen abundances from the IR OH lines are obviously too high 
with respect to what would be expected for the EMP stars in this 
[Fe/H] rang43- On the other hand, 3D LTE abundances obtained 
from OH lines agree significantly better with the measurements 
obtained from [Oi] line than the ID LTE estimates do (Fig.|4l). 


We note that similar trends of oxygen abundance with 
equivalent width w e re als o no ticed in the i nvestig ations of 
iMelendez & Barblr^ (l2002l) and iBarbuv et alJ (l2003l) and they 
remain unexplained. We also investigated if the determined oxy¬ 
gen abundances depend on quantum number J, but could find 
no clear correlation. It is encouraging that, ignoring the trends, 
the line-to-line scatter of the abundances derived from the differ¬ 
ent lines is small enough so that the mean can be considered a 
reasonable estimate of the oxygen abundance in a given star. 

One can conclude that there are hints of a slight dependence 
of abundance on excitation potential in the case of BD-18:5550 
(and perhaps also CS 31081-001; see Fig. |2]i. Nevertheless, in 
the case of BD-18:5550 the two lines with the highest excitation 
potential and the highest oxygen abundance also significantly 
deviate in the abundance-equivalent width plane (Fig. [TJ. Why 
abundances obtained from these two lines differ is not clear; we 
note, however, that these lines are the weakest ones measured in 
the spectrum of this star. 

We identify two ways to improve our modelling of these 
lines. The first is to increase the number of 3D model atmo¬ 
spheres used in the analysis. One clear limitation of the present 
investigation is the use of just one model atmosphere to describe 
the four programme stars, that, although they have been selected 
with very similar atmospheric parameters, span a range in sur¬ 
face gravities, effective temperatures, and metallicities. In par¬ 
ticular, the sensitivity of the 3D corrections to the variations of 
atmospheric parameters needs to be investigated. Such an inves¬ 
tigation is hampered by the long time needed to fully relax a 3D 
model atmosphere of a giant star. 

The second way is the investigation of possible departures 
from the assumption of LTE, that has been made in the present 
investigation. Departures could occur both in the chemical equi¬ 
librium (formation and dissociation of OH molecules here are 
assumed to be an equilibrium process); and in the occupation 
number of different energy levels of the OH molecule (the Boltz¬ 
mann law was assumed here). 

In Fig.|4]we show the [O/Fe] ratios in the four programme 
stars together with data taken from the literature based on the 
630 nm [Oi] line. The oxygen abundances are the ones we 
determined from the IR OH line s and the iron abundances 
have been taken from ICavrel et al.l (l2004ll . The points labelled 
3D mean that only the oxygen abundance has been corrected 
for the 3D effect, but not the iron. In the case of iron, the 
3D correc tion for such a mode l atmosphere should be ar ound 
-0.2 dex (iBonifacio et alJ 120091: iDobrovolskas et ak 2013h. al¬ 
though strongly dependent on excitation potentia l ( Collet et~aD 
l2007tl^llet et alJlioOOtfDobrovolskas et alJl2013h . On the other 
hand, we also know that Fei lines imply NLTE corrections 
of the same order of magnitude but in the opposite direc¬ 
tion. In thei r comprehensive s t udy o f NLTE effects on iron in 
HD 122563 iMashonkina et akl (|2()1 lb determined a NLTE cor¬ 
rection of - 1 - 0 .16 dex (for collisions with neutral hydrogen treated 


The [O/Fe] abundance ratios from ICavrel et all l l2004h were obtained 
using solar oxygen abunda nce of A(0) = 8.74, whereas in our case 
we used A(0) = 8.66 from lAsplund et aP (|2005h . To m ake the direct 
comparison possible , the [0/Fe1 rat i os fro m ICavrel et all (l2004h and the 
theoretical model of iFrancois et all ( 120041) were shifted to the reference 
scale where the solar abundance is A(0) = 8.66. 


with the Drawin formalism and an efficiency Sh = 0.1; neglect¬ 
ing the coll isions raises the NLTE co rrection to +0.35). On the 
other hand, IMashonkina et al.l (12013l) demonstrated that the 3D 
and NLTE corrections for iron cannot simply be added; this ap¬ 
proach for HD 122563 leads to a discrepancy of -0.75 dex be¬ 
tween the abundance derived from the Fei and Feii lines. The 
abundance of iron probably requires a full 3D-NLTE analysis, 
and we are currently unable to perform it. In this state of affairs 
we prefer to adopt the ID LTE result for iron. 

With this choice, the 3D [O/Fe] values obtained from OH 
lines appear to be in better agreement with the general trend de¬ 
fined by the abundances based on [Oi]. The average 3D LTE 
oxygen abundance for the four stars obtained from OH lines is 
{A(0)) = 6.61 + 0.21, while that obtained from the 630 nm line 
is (A(0)) = 6.52 + 0.33. We note that in ID LTE we obtained 
{A(0)) = 6.86 ± 0.22 from OH lines; errors in all cases here 
indicate the scatter due to star-to-star abundance variations. 

The results shown in Fig. |4] are, by and large, con- 
sistent with those deriv e d fro m the OH UV lines by 
[Gonzalez Hernandez et al.l (1201 Oh in metal poor dwarfs, al¬ 
though these authors have adopted the strategy of adding 3D 
corrections and ID NLTE correction for Fei. This means that 
a direct comparison of the two results is not possibl e. Never¬ 
theless, the results of I Gonzalez Hernandez et alJ (1201 Ol see their 
Fig. 14) show an average [O/Fe] = 0.6 - 0.8, corrected for 
3D and NLTE effects at [Fe/H] = -3.0. In our case, we ob¬ 
tain an average [O/Fe] = 0.7, corrected only for 3D effects, at 
[Fe/H] = -3.0. On the other hand, the r esults of the two stud- 
ies agree well with theoretical model of iFrancois et al.l (l2004h 
which predicts [O/Fe] = 0.8 for the same metallicity (in this 


co mparison, [O/Fe 


value s from iGonzalez Hernandez et alJ2010| 


and [ Francois et al]|2004l were recomputed to the lAsplund et al. 


12005 solar oxygen abundance scale). In the future we plan to 
investigate the consistency of abundances obtained from IR OH 
and UV lines in the same EMP stars. 


5. Conclusions 

The results obtained in our pilot project have demonstrated that 
infrared OH vibrational-rotational lines could be used success¬ 
fully to determine oxygen abundance in EMP giants. This, how¬ 
ever, requires high-resolution spectra obtained with S IN ratios of 
at least a few hundred. Fortunately, IR OH lines are significantly 
stronger than the [O i] 630 nm line and thus can be more easily 
measured. In addition, quite a large number of IR OH lines un¬ 
affected by telluric lines are typically available for the measure¬ 
ments, which makes the statistical significance of the determined 
oxygen abundance higher. 

To determine oxygen abundances from IR OH lines in EMP 
stars, however, the use of 3D hydrodynamical model atmo¬ 
spheres is imperative because these lines form in the outer at¬ 
mosphere where horizontal temperature fluctuations (caused by 
convective motions) are large, which leads to large variations in 
line opacities. This results in significantly different line strengths 
predicted with the 3D hydrodynamical and ID hydrostatic model 
atmospheres, with the resulting 3D-1D differences in oxygen 
abundances of -0.2 ■ ■ ■ - 0.3 dex. When corrected for these ef¬ 
fects with the aid of 3D hydrodynamical model atmospheres, 
oxygen abundances obtained from IR OH lines generally agree 
well with those determined using the [O i] 630 nm line, which is 
insensitive to either 3D or NLTE effects. Obviously, the role of 
NLTE effects in the IR OH line formation still needs to be clar¬ 
ified, as well as the combined 3D NLTE effects for iron, some¬ 
thing that is anticipated from the forthcoming studies. 
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Table B.l. Abundance corrections for HD 122563. 


Table B.3. Same as in Table lB.ll but for BD-18:5550. 


nm 

EW, 

(pm) 

3D - <3D) 
(dex) 

(3D) - ID 
(dex) 

3D-1D 

(dex) 

d, 

nm 

EW, 

(pm) 

3D - (3D) 
(dex) 

(3D) - ID 
(dex) 

3D-1D 

(dex) 

[O i] line 





[Oi] line 





630.0.304 

0.64" 

-0.044 

0.039 

-0.005 

630.0.304 

0.15" 

-0.043 

0.038 

-0.006 

OH lines 





OH lines 





1514.7939 

1.20 

-0.41 

0.13 

-0.29 

1514.5770 

0.515 

-0.41 

0.12 

-0.29 

1527.8524 

1.30 

-0.41 

0.12 

-0.28 

1514.7939 

0.562 

-0.41 

0.12 

-0.29 

1528.1052 

1.38 

-0.41 

0.12 

-0.28 

1527.8524 

0.560 

-0.41 

0.12 

-0.28 

1605.5464 

1.09 

-0.34 

0.12 

-0.21 

1528.1052 

0.636 

-0.41 

0.12 

-0.28 

1619.0132 

1.71 

-0.33 

0.12 

-0.21 

1606.1700 

0.390 

-0.36 

0.12 

-0.24 

1619.2130 

0.97 

-0.33 

0.12 

-0.21 

1606.9524 

0.361 

-0.36 

0.12 

-0.24 

1620.4076 

1.84 

-0.33 

0.12 

-0.21 

1623.0465 

0.352 

-0.30 

0.12 

-0.18 

1625.5019 

1.61 

-0.35 

0.12 

-0.23 

1623.1419 

0.349 

-0.30 

0.12 

-0.18 

1626.0155 

1.63 

-0.35 

0.12 

-0.23 











^ adopted fromICavrel 

^tlildTOOl 



^ adoDted fromICavrel 

etal. 1(12004 













Table B.4. Same as in TablelObut for CS 31082-001. 


fable B.2. Same as in Table|BT]but for HD 186478. 












4, 

EW, 

3D - (3D) 

(3D) - ID 

3D-1D 

4, 

EW, 

3D - <3D) 

(3D) - ID 

3D-1D 

nm 

(pm) 

(dex) 

(dex) 

(dex) 

nm 

(pm) 

(dex) 

(dex) 

(dex) 

[Oi] line 





[O i] line 





630.0.304 

0.27" 

-0.044 

0.038 

-0.006 

630.0.304 

0.93" 

-0.045 

0.040 

-0.005 

OH lines 





OH lines 





1526.6167 

0.82 

-0.41 

0.12 

-0.28 

1514.5770 

1.52 

-0.42 

0.13 

-0.29 

1527.8524 

0.89 

-0.41 

0.12 

-0.28 

1514.7939 

1.71 

-0.42 

0.12 

-0.29 

1619.0132 

0.97 

-0.33 

0.12 

-0.21 

1526.4604 

1.13 

-0.41 

0.12 

-0.28 

1620.4076 

0.77 

-0.33 

0.12 

-0.21 

1526.6167 

1.46 

-0.40 

0.12 

-0.28 






1527.8524 

1.60 

-0.41 

0.12 

-0.28 

^ adopted fromlCayrel 

^tlfldTOOl 



1528.1052 

1.53 

-0.41 

0.12 

-0.28 






1541.9460 

1.69 

-0.40 

0.12 

-0.27 






1542.2366 

1.87 

-0.39 

0.12 

-0.27 






1605.2765 

1.44 

-0.34 

0.12 

-0.22 






1606.1700 

1.39 

-0.36 

0.12 

-0.24 






1625.1660 

0.93 

-0.35 

0.12 

-0.23 






1625.5019 

1.61 

-0.35 

0.12 

-0.23 






1626.0155 

1.22 

-0.35 

0.12 

-0.23 







“ adopted from ICavrel et al.l (l2004h 


Appendix A: List of infrared OH iines for oxygen 
abundance determination 

In this work we used a number of vibrational-rotational lines 
(A^^n) from the first-overtone sequence (An = 2) located in the 
spectral range between 1514 - 1548 and 1595 - 1632 nm. The 
atomic parameters of the IR OH lines and their equivalent widths 
measured in the target EMP giants and also the derived oxygen 
abundances are provided in Table I A. fl 


Appendix B: 3D-1D abundance corrections for 
iR OH iines and the forbidden oxygen iine 

In Tables lBTTIIB.4l we provide 3D-1D abundance corrections that 
were computed for all OH lines and forbidden oxygen line, in 
each EMP giant studied here. 
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Table A.l. The list of the OH spectral lines used in the abundance analysis. 


nm 

A. 

eV 

log 9/ 

v' - v" 

HD 122563 
IV, A(0) 
(pm) 

HD 186478 
IV, A(0) 
(pm) 

BD-18;5550 
IV, A(0) 
(pm) 

CS 31082-001 
W, A(0) 

(pm) 

1514.5770 

0.164 


0-2 



1.52 

7.15 

0.52 

6.72 



1514.7939 

0.164 

-5.447 

0-2 

1.20 

6.61 

1.71 

7.20 

0.56 

6.76 



1526.4604 

0.210 

-5.428 

0-2 



1.13 

7.04 





1526.6167 

0.210 

-5.428 

0-2 



1.46 

7.16 



0.82 

7.09 

1527.8524 

0.205 

-5.382 

0-2 

1.30 

6.62 

1.60 

7.15 

0.56 

6.72 

0.89 

7.07 

1528.1052 

0.205 

-5.382 

0-2 

1.38 

6.64 

1.53 

7.12 

0.64 

6.78 



1541.9460 

0.250 

-5.323 

0-2 



1.69 

7.15 





1542.2366 

0.250 

-5.323 

0-2 



1.87 

7.06 





1605.2765 

0.639 

-4.910 

1-3 



1.44 

7.09 





1605.5464 

0.640 

-4.910 

1-3 

1.09 

6.54 







1606.1700 

0.476 

-5.159 

0-2 



1.39 

7.11 

0.39 

6.60 



1606.9524 

0.472 

-5.128 

0-2 





0.36 

6.53 



1619.0132 

0.688 

-4.893 

1-3 

1.71 

6.57 





0.97 

6.91 

1619.2130 

0.688 

-4.893 

1-3 

0.97 

6.52 







1620.4076 

0.683 

-4.851 

1-3 

1.84 

6.69 





0.77 

6.84 

1623.0465 

0.905 

-5.059 

2-4 





0.35 

6.81 



1623.1419 

0.905 

-5.059 

2-4 





0.35 

6.90 



1625.1660 

0.543 

-5.115 

0-2 



0.93 

6.96 





1625.5019 

0.538 

-5.087 

0-2 

1.61 

6.75 

1.61 

7.17 





1626.0155 

0.539 

-5.087 

0-2 

1.63 

6.75 

1.22 

7.04 














normalized flux 
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HD 122563 





— •-CRIRES - Gaussian fit - TAP AS 








Fig. B.l. Excerpts from CRIRES spectra centred on IR OH lines in HD 122563 (black dots). Gaussian best-fits to the o bserved IR OH line p rofiles 
are shown as red lines. Blue solid lines show the synthetic telluric spectrum computed with the TAPAS online service dBertaux et al.ll2014li which 
is shifted downwards by 0.01 for better readability. Rest frame wavelengths (do) and measured equivalent width values {EW) of the IR OH lines 
are provided at the bottom of each panel. 
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HD 186478 





— • —CRIRES - Gaussian fit - TAP AS 


T- 1 - 1 - 1 -r 






1514.7942 nm 
EW= 1.71 pm 

_l_I _ I _I_ I _ I_L 


1514.8 


1515.0 







Fig. B.2. The same as in Fie. lB.ll but for HD 186478. 


























normalized flux 


HD 186478 
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— • — CRIRES - Gaussian fit - TAP AS 



, I - 1 - 1 - 1 - 1 - 1 - 1 - 


1625.1662 run 
- £'lF=0.93pm 

^^^ _ 


1625.0 


1625.2 


1 - 1 -r 


1.00 






0.95 - 

1626.0154 nm 

- £»"=1.22pm 

0.90 -1-1-1-1-1-1-L 


1625.8 


1626.0 

T, nm 


1625.4 



1626.2 


Fig. B.3. The same as in Fie. IB. 21 












BD-18:5550 
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—CRIRES - Gaussian fit - TAP AS 



\ = 1514.7942 run 
£'lF=0.56pm 


± 


J_I_L. 


± 


1514.0 


1514.2 





1.00 


0.90 


- ^ ^ ^ ^ ^ ^ ^ - 


1 1 1 1 1 1 1 

- 2^ = 1623.0455 nm 


2^ = 1623.1415 nm 

- FlF=0.35pm 


. FfT=0.35pm 

_1_1_1_1_1_1_1_ 


_1_1_1_1_1_1_1_ 


1622.2 


1622.4 


1622.2 


1622.4 
/I, nm 


1622.6 


Fig. B.4. The same as in Fie. lB.ll but for BD -18:5550. 
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CS 31082-001 


— • — CRIRES - Gaussian fit - TAP AS 


0.95 - 


0.90 


-1-1-1--1-1-1- 


]T ' ^ ^ ^ " 


-1-1-1-1-1-1-1- 

— JFV y 





] 1526.6167 nm “ 


1527.8524 nm “ 


- /i^= 1619.0132 nm 

FW- 0 Q7 r>m 

EW- 0.82 pm 

_1_1_1_1_1_1_1_ 


EW= 0.89 pm 

_1_1_1_1_1_1_1_ 


_1_1_1_1_1_1_1_ 


1527.0 


1527.2 


1527.4 


1528.4 


1528.6 


1619.6 


1619.8 


1.00 


0.95 - 


0.90 


1-1-1-1-r 




= 1620.4076 nm 
£'r=0.77pm 

_J_I _ I _I_ L_ 


1620.8 1621.0 

/I, nm 


1621.2 


Fig. B.5. The same as in Fie. lB.ll but for CS 31082-001. 














